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Abstract. In this paper we systematically apply the mathematical structures by time-evolving 
measures developed in a previous work to the macroscopic modeling of pedestrian flows. We 
propose a discrete-time Eulerian model, in which the space occupancy by pedestrians is described 
via a sequence of Radon positive measures generated by a push-forward recursive relation. We 
assume that two fundamental aspects of pedestrian behavior rule the dynamics of the system: 
On the one hand, the will to reach specific targets, which determines the main direction of motion 
of the walkers; on the other hand, the tendency to avoid crowding, which introduces interactions 
among the individuals. The resulting model is able to reproduce several experimental evidences 
of pedestrian flows pointed out in the specialized literature, being at the same time much easier 
to handle, from both the analytical and the numerical point of view, than other models relying 
on nonlinear hyperbolic conservation laws. This makes it suitable to address two-dimensional 
applications of practical interest, chiefly the motion of pedestrians in complex domains scattered 
with obstacles. 



1. Introduction 

In recent years, multi-agent systems, like e.g., cell populations, bird and fish swarms, insect 
colonies, human crowds, have roused the interest of experts in various fields, from biology to 
sociology, physics, and applied mathematics. The primary reason for this scientific attention 
has probably to be sought in that such systems issue new stimulating challenges toward their 
comprehension. Indeed, they are strongly nonstandard, because of some forms of intelligence and 
decisional abilities, which can affect in a crucial way their evolution. In addition, they are complex 
systems, in which group dynamics plays an essential role: The collective behavior is not the simple 
superposition of individual behaviors, for complex interactions usually arise among the subjects, 
leading finally to completely new dynamics with respect to the case of isolated agents (cf. Krause 
and Ruxton [24]). 

Design and use of innovative analytical and numerical tools to address the study of these systems 
is one of the big challenges of modern applied mathematics. This is motivated, first of all, by the 
interest in stating a qualitative and quantitative formalization of the above-mentioned behaviors, 
which may help improve traditional socio-biological investigation methods by taking advantage 
of the typical conciseness of mathematics. On the other hand, referring in particular to human 
crowds, numerous engineering applications welcome the support of mathematical modeling in this 
nonclassical field for design and optimization purposes. For instance, in the structural engineering 
circle it is well-known the case of the London Millennium Footbridge, closed the very day of 
its opening due to macroscopic lateral oscillations of the structure developing while pedestrians 
crossed the bridge (cf. Dallard et al. [9]). This unexpected phenomenon, caused by a synchronous 
lateral excitation exerted by the walkers on the footbridge, pointed out the possible relevance of 
crowd- structure coupling in the design of walkway infrastructures, and renewed the interest toward 
the investigation of these issues by means of mathematical modeling techniques (cf. Venuti et al. 
[30] and the main references therein). Another example is the Jamarat Bridge, located in Mina, 
South Arabia, 5 km far from Mecca, reached each year by a huge crowd of pilgrims for the Hajj, 
the pilgrimage to Mecca that, according to the prescriptions of Islam, Muslims must carry out at 
least once in their lifetime. The extraordinary amount of people cramming the bridge in those 
occasions gave rise to serious pedestrians disasters in the nineties, with a lot of people died for 
overcompression. In this case, mathematical modeling and numerical simulation have already 
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been successfully employed to study the dynamics of the flow of pilgrims, so as to highlight critical 
circumstances under which crowd accidents tend to occur and suggest counter-measures to improve 
the safety of the event (cf. Helbing et al. [12], Hughes [22]). 

The experimental literature on pedestrian behavior is already quite evolved. Observations and 
data recording by the most sophisticated techniques are nowadays extensively available [lOj [TTJ 
[12] [15] [19] [20] . thanks to which numerous human walking attitudes have been pointed out. For 
instance, it is known that pedestrians tend to maintain a preferential direction of motion toward 
some targets they want to reach, but at the same time they are disposed to slightly deviate from it 
in order to avoid crowding. In evaluating the amount of people occupying the neighboring space, 
walkers take into account almost exclusively what they see ahead, because the particular position 
of the eyes on the head reduces their visibility field to a frontal area only. This has important 
consequences on the spatial configurations assumed by moving pedestrian masses: Pedestrians 
tend to line up in parallel lanes, unlike other kinds of agents, e.g.. birds and fishes, which instead 
prefer to flock. Notice that the visibility field of the latter covers also a rear area, due to the lateral 
positioning of their eyes on the head. Even more interesting is the well-studied case of oppositely 
walking pedestrian groups, in which the above-mentioned lane formation is particularly evident: 
When the groups cross each other, the emergence of alternate uniformly walking lanes is observed 
as a characteristic self-organizing tendency. 

Like for many other real-world systems, collecting huge amounts of experimental data can be 
only a first step toward the scientific investigation of pedestrian flows. Indeed, phenomenological 
observations certainly enable one to gain fundamental insights into typical aspects of this complex 
system, as exemplified above, but they hardly allow to catch the dynamics far from equilibrium 
conditions. Consequently, they are mainly descriptive but scarcely predictive, as the unsteady 
dynamics often plays a relevant role in determining the overall evolution. Crowd motion is therefore 
definitely not only an engineering matter but also a challenging mathematical problem, which calls 
for the development of efficient modeling and simulation tools. Unlike classical physical systems 
dealing with inert matter, a well-coded mathematical- socio-biological theory is so far still lacking, 
whence the current impossibility to derive mathematical models of the behavior of the multi-agent 
systems discussed in this introduction from first principles. Therefore, specific models will feature 
necessarily phenomenological aspects. However, mathematical structures can be developed in the 
abstract by appealing to the most sophisticated techniques, with the aim of providing efficient 
modeling frameworks from both the theoretical and the applied point of view. 

In this paper we are concerned with mathematical modeling of pedestrian flows. Five more 
sections follow this introduction, namely: 

• Section [2] in which we account for the main literature on mathematical models of pedes- 
trian flows currently available. Our short review covers models at both the microscopic 
and the macroscopic scale. 

• Section [3] in which we illustrate the mathematical structures by time-evolving measures 
that we have introduced in our previous work [28 for the macroscopic modeling and the 
numerical simulation of continuous systems. We also summarize some related theoretical 
results, however without insisting on their proofs. The interested reader can find all details 
in our above-cited paper. 

• Section[4] in which we specifically apply the above framework to the modeling of pedestrian 
flows, discussing in particular the structure of the velocity field of pedestrians and the 
handling of boundary conditions. The latter issue pertains in particular to the modeling 
of the interactions between pedestrians and walls or obstacles. 

• Section [5] in which we address some relevant cases study by the model previously derived. 
We consider the motion of pedestrians in areas scattered with obstacles by simulating a 
group of travelers accessing one or two escalators after getting off a subway, and a group 
of walkers wanting to get out of a room in which two pillars partially hide the exit. Then 
we focus on self-organization phenomena, showing that our model is able to reproduce 
the typical lane formation and the emergence of alternate lanes in crossing flows. By 
slightly modifying a modeling detail, so as to account for a full ahead-rear visibility area 
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of the agents, we also suggest that our model is able to qualitatively reproduce flocking 
phenomena, but we refrain from going too much into this issue as it is beyond the scope 
of the present work. 

• Section [6j in which we draw some conclusions and briefly sketch research perspectives. 

2. Overview of mathematical models of human crowds 

In this section we briefly account for the existing literature on mathematical modeling of pedes- 
trian flows, in order to depict the scientific frame which the present paper fits in. Despite the 
relative novelty of the subject in the area of applied mathematics (first microscopic and macro- 
scopic models due to Helbing, Hughes, and their respective coworkers date back to the late nineties 
and early two-thousands), several relevant contributions to this research line can be found already, 
including a number of works dealing with experimental investigation and data recording on human 
walking attitudes. 

Pedestrians in motion within a given walking area are an essentially discrete system, in which 
each person plays the role of an isolated agent moving in interaction with other surrounding agents, 
usually with the aim of reaching a particular destination. As such, this system can be conceptually 
described at the microscopic scale, i.e., by a system of ordinary differential equations tracing the 
evolution in time of the spatial position of each single pedestrian. The possible coupling of these 
ODEs should reflect the interactions among pedestrians, namely the influence that each one of 
them has on the motion of the others. Some microscopic models of pedestrian flow are currently 
available in the literature, generally regarding pedestrians as rigid particles, whose motion is 
regulated by classical Newtonian laws of point mechanics: 

Vi(t)=Fi(t), 

where V{ is the velocity of the i-th pedestrian, Fi an overall force acting on her/him, while the 
index i runs from 1 to the total number N of pedestrian considered in the model. The analogy 
with standard mechanical systems of inert matter is, however, only formal, as the force fields Fi 
invoke nonclassical dynamic concepts resorting mainly to the ideas of: 

• Preferred direction of motion, i.e., the fact that pedestrians normally have a target to walk 
toward, which guides the main direction of their movement. 

• Comfort/ discomfort at certain distances from other pedestrians, i.e., the fact that in nor- 
mal conditions pedestrians may feel uncomfortable if too close to one another and tend to 
look for uncrowded surrounding areas. 

Therefore, in this context the F^s are not understood strictly as mechanical actions exerted by 
pedestrians on each other. For instance, Helbing and coworkers [T3l E] introduce the concept of 
social (or behavioral) force, which measures the internal motivation of the individuals in performing 
certain movements. Specifically, in their model pedestrians are regarded as points, and two main 
factors contribute to the definition of the social force Fi acting on the z-th individual: 

• A relaxation toward a desired velocity v®, i.e., the velocity that the i-th individual should 
possess in order to reach her/his destination as comfortably as possible: 



Ti > being the characteristic relaxation time of pedestrian i. 
• A repulsive effect from either neighboring pedestrians located too close, i.e., within the 
so-called private sphere of the i-th pedestrian: 

fij( r ij) = ~^ V ij( r ij)^ 3 = 1> • • • j N , 

where Vij is a monotonically decreasing repulsive potential from pedestrian j with ellipse- 
shaped level curves turned in the direction r^ = Xj — Xi, or from edges of walls and 
obstacles found within the walking area: 

F iE (r iE ) = -VU iE (\\r iE \\), 
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where Uie is a monotonically decreasing repulsive potential from edge E and = x l E — 
x l E G E being the edge's nearest point to pedestrian i. 

After setting Fi = Ff + YljLi fij + ^2e FiE, where the last sum at the right-hand side is extended 
to all edges E of walls and obstacles that may affect the motion of the individuals, the final system 
of equations reads: 



dvj _ - Vj _ y> 
dt ~ n ^ 



V^(r^) - VU iE (\\r iE \\), i = 1, . . . , N. 



A more refined version of the model includes, still using suitable scalar potentials, also possible 
attractive effects exerted on pedestrians by particular elements such as windows, displays, other 
people, as well as stochastic fluctuations in the behavior of each subject. The interested reader is 
referred to Helbing and Molnar [15] , and to the main references listed therein, for further details 
on these issues. 

Similar behavioral considerations underlie the microscopic model of pedestrian flow by Maury 
and Venel [27], although its formalization resorts to different mechanical guidelines. Pedestrians 
are still regarded as rigid particles, in particular disks of fixed radii, but the emphasis is now on 
a kinematic construction of their actual velocity, which does not invoke explicitly any concept of 
generalized force. In more detail, the preferred direction of motion is expressed by introducing 
a desired velocity J7, determined by the geometry of the walking area in such a way that it 
drives pedestrians toward their target along the shortest path, taking into account at the same 
time the possible presence of intermediate obstacles. On the other hand, the repulsive effect 
among the individuals is meant as a geometrical constraint to avoid that disks step over one 
another. Specifically, it is fulfilled by first defining a set C q of admissible velocities, which at 
each time instant depends on the positions, denoted by ^, of all pedestrians i = 1, . . . , iV, and 
then constructing the actual velocity as the Euclidean projection of the desired velocity of each 
pedestrian onto C q . The evolution in time of the system is finally described as 



Qi(t) = <fc |0 + J (P Cq U)(s)ds, i = 1, 



o 

where q iy o is the initial position of pedestrian i and Pc q is the Euclidean projector onto C q . 

A further microscopic formulation of the problem is instead used by Hoogendoorn and Bovy 
[16] [18] , who propose a theory of pedestrian behavior based on the concepts of walking task 
and walking cost. Basically, they assume that pedestrians are feedback-oriented controllers, who 
plan their movements on the basis of some predictions they make on the behavior of the other 
individuals. Predictions are dictated by a sort of cooperative or non-cooperative game theory. 
In either case, they are affected by a limited in time and space predictive ability of the walkers. 
Each pedestrian behaves so as to minimize her/his individually estimated walking cost, which is 
expressed by a suitable functional depending on the predicted positions of the other people. 

A different approach to the description of the system uses instead partial differential equations 
and the theory of (possibly multidimensional) conservation laws. In this case, it is assumed that 
pedestrians moving within a given walking area have a continuous distribution in space, so that it 
makes sense to introduce their density p = p(t, x) and to invoke some conservation principles, e.g. 
the conservation of mass and possibly also of linear momentum, in order to get an equation, or a 
system of equations, satisfied by p. One of the first threads of macroscopic models of pedestrian 
flow is due to Hughes [2TJ [22] [23] , who proposes a two-dimensional diffusion-like equation 

- V • ( W (p)/ 2 (p)W) = (2.1) 

coming from the mass conservation equation supplemented by the following phenomenological 
closure of the velocity: 

v = -g{p)f{ P )Vct>. 

Here, <fi is a scalar potential whose gradient V0 determines geometrically, at each point of the 
domain, the main direction of motion of the individuals, whereas /, g are scalar factors affecting 
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the speed of pedestrians on the basis of the local crowding of the walking area, in particular 
considering that pedestrians try to avoid extremely high densities. By means of model (2.1) the 
Author studies various two-dimensional configurations of human crowd motion, including the flow 
of pedestrians past obstacles with application to the celebrated case of Jamarat bridge (see also 
Helbing et at. [12 for an experimental investigation of this issue). 

Colombo and Rosini [6i introduce instead a one-dimensional macroscopic model built on a 
Cauchy problem for the nonlinear hyperbolic conservation law 

d t p + d x f(p) = 0, xGl,t>0, 

which reminds of the Lighthill-Whitham-Richards (LWR) model of vehicular traffic (see Lighthill 
and Whitham [25], Richards [29]). The main difference is that the density of pedestrians exhibits 
two characteristic maximal values R, R* , with < R < R*, at both of which the flux / vanishes. 
In normal situations p ranges in the interval [0, R], where the flux is nonnegative, either strictly 
concave or with at most one inflection point, and has precisely one local maximum point. When 
the density grows above its "standard" maximum value, i.e., for p G (i?, i?*], it is assumed that, 
unlike vehicular traffic, pedestrians can still move but feel overcompressed, hence their flux is 
less effective than before and they enter a panic state. In this region, the function / features 
a trend similar to that described for p G [0, R], but with a local maximum value strictly less 
than the previous one. As illustrated in [6j[7], this allows to define a concept of solution to the 
above conservation law in which non-classical shocks are admitted, i.e., shocks complying with 
the Rankine-Hugoniot condition but possibly violating entropy criteria. As a consequence, the 
classical maximum principle for nonlinear hyperbolic equations, stating that the solution p(t, x) 
remains confined within the same lower and upper bounds of the initial datum for all x G M and 
all t > 0, no longer holds true and the model is able to describe the transition of pedestrians 
to panic even starting from an initial density entirely bounded below the standard maximum R. 
The resulting fundamental diagram, i.e., the mapping p /(p), agrees well with experimental 
observations reported by Helbing et al. in [12] . 

Bellomo and Dogbe [3], Coscia and Canavesio [8] refer instead to a two-dimensional setting, 
in which the walking area is represented by a bounded domain Q C M? with possible inlet and 
outlet regions along the boundary dfl. In [3 the motion of pedestrians is described by a system 
of two partial differential equations invoking the conservation of mass and the balance of linear 
momentum: 

f^p + v-(H = o 

\dtv + (vV)v = F\p,v], [ ' } 

where F is a material model for the acceleration of the individuals depending in general in a 
functional way on the density p and the velocity v. Conversely, in [8] only the continuity equation 
(first equation in (2.2)) is used, and the model is made self-consistent by devising appropriate 
closure relations for the velocity v in terms of the density p and possibly also of its gradient Vp: 

v\p, Vp](x) = <p\p, Vp]i/(x), 

where v = v{x) is a unit vector identifying, at each point x G the preferred direction of motion 
of the crowd and ip a scalar function expressing the speed of the latter (square brackets denote in 
this context functional dependence). The Authors suggest several forms for <p, referring sometimes 
also to the velocity diagrams of first-order vehicular traffic models. 

Equations (2.2) are formally inspired by the classical fluid dynamics models of continuum 
mechanics, however the force (per unit mass) F contains non-classical contributions accounting 
for: 

• A relaxation toward a desired velocity, that makes pedestrians point in the direction of a 
certain target they want to reach: 

F x = a(v e (p)v - v), 

where a > is the inverse of the relaxation time of pedestrians, vq is a unit vector pointing 
toward the target, and finally v e (p) is an equilibrium speed of pedestrians depending 
pointwise on the crowding of the walking area. 
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• A local crowding estimate, based on the pointwise values that Vp takes along the direction 
of the desired velocity, which might induce pedestrians to deviate from their preferred path 
in order to avoid areas of high density: 

^ = K\p) dp ^ 
P du ' 

A pressure-like term, possibly regarded as a material quantity as in the celebrated Aw- 
Rascle model of vehicular traffic pQ, which models the reaction of pedestrians to the 
presence of other individuals in the surrounding environment. When also this term is 



included, the linear momentum balance (second equation in (2.2)) rewrites as 



d t (v + P(p, v)v ) + (v ■ V)(v + P(p, v)v ) = F[p, v], 

where P = P(p, v) is some pressure that walkers feel along the preferred path, depending 
on the pointwise crowding of the domain and on their current velocity. 

Additional topics, like e.g. the existence of a limited visibility zone for each pedestrian when trying 
to evaluate the minimal crowding direction, are also discussed. In particular, special attention is 
paid to the characterization of the panic state and to the transition to it from regular conditions: 
The Authors suggest that pedestrians entering a panic state tend to follow chaotically other 
individuals, dropping any specific target, and therefore are mostly attracted toward areas of high 
density rather than seeking the less congested paths. 

The usually large amount of coupled ODEs to be handled simultaneously, that microscopic 
models usually require, is a drawback if not from the computational point of view, thanks to the 
increasing power of modern calculators, certainly for analytical purposes. Among others, we recall 
here the difficulty to recover a global overview of the system from the knowledge of its microscopic 
state, possibly in connection with control and optimization issues. Macroscopic models are more 
suited to this, but those currently available in the literature have to face several other complications 
due to their intrinsic hyperbolic nature. First of all, the most natural settings for pedestrian flow 
problems are two-dimensional: One-dimensional models are essentially explorative, but they are 
unlikely to provide effective mathematical tools to deal with real applications. However, it is 
well known that the theory of multi-dimensional systems of nonlinear hyperbolic equations is 
much more complicated under both the analytical and the numerical point of view, therefore a 
sound mathematical mastery of such models may hardly be achieved. Secondly, the imposition 
of boundary conditions may be tricky in hyperbolic models, because on the one hand one is 
forcedly driven by the characteristic velocities in defining the inflow and the outflow portion of 
the boundary, while on the other hand it must be guaranteed that pedestrians do not enter or exit 
the domain from any point of the boundary other than the prescribed inlet and outlet regions. 
This issue gets even more complicated in presence of obstacles, which have to be understood as 
internal boundaries to the walking area. We notice, however, that the most interesting problems 
for applications generally do not concern pedestrian motion in free spaces, but precisely in areas 
scattered with obstacles (e.g., pillars, bottlenecks, narrow passages, see Helbing et al. fM]). 
sometimes used to force the flow of crowds in specific directions. 

A time-evolving measures approach to the modeling of pedestrian flow may help overcome some 
of the difficulties just outlined. As we will see in much more detail in the next section, the basic idea 
is to use measures as mathematical tools to evaluate the degree of space occupancy by pedestrians. 
In this way, it is perfectly natural to address the problem from a macroscopic, even Eulerian, point 
of view in spite of the intrinsic Lagrangian granularity of the system: Given a walking area 
one describes the evolution of the system by measuring the crowding of every subset E C £1 at 
each time instant. In addition, as we have shown in [28] , there are basically no differences in the 
one- or two- (and even three-) dimensional theory, therefore one can immediately tackle realistic 
problems without the need for conceiving preliminary one-dimensional approximations. 
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3. Mathematical modeling by time-evolving measures 

Canuto and coworkers have proposed in [5] an Eulerian measure-theoretical approach to the 
modeling of rendez-vous problems for multi-agent systems. Inspired by their work, we have intro- 
duced in [28 a measure-theoretical continuous modeling framework based on conservation laws, 
which then has been preliminarily applied to pedestrian flow problems. The main novelty with 
respect to the standard literature on the subject is that such conservation laws are stated in terms 
of discrete-time-evolving measures rather than through hyperbolic partial differential equations, 
which yields several advantages over classical approaches also used by some Authors to face similar 
problems. From the analytical point of view, we recall: (i) An easier passage from the Lagrangian 
to the Eulerian description of the system, without the need for resorting heavily to regularity 
issues; (ii) A more straightforward establishment of the well posedness of the problem in terms of 
existence, uniqueness, and a priori estimates of the solution; (iii) A direct deduction of an ad-hoc 
numerical scheme for the approximate treatment of the equations, characterized by nice stability 
and accuracy properties as well as by a simple practical implementation. From the modeling point 
of view, we anticipate that this framework: (iv) Allows to address directly realistic two-dimensional 
applications without extra difficulties with respect to the (explorative) one-dimensional case; (v) 
Is particularly suited to treat nonlocal interactions among pedestrians; (vi) Gives rise to an easy 
and rich handling of boundary conditions as far as both wall-like and obstacle-like boundaries are 
concerned. 

3.1. Modeling framework. Let ft C M d , where d — 1, 2, 3 from the physical point of view, 
be a bounded set representing the walking area of pedestrians. Topologically, we may think of 
it as a pathwise-connected domain possibly containing holes understood as internal obstacles to 
the walking area. The core of the modeling approach is the description of the space occupancy 
by pedestrians via a family of positive Radon measures {/i n }n>o defined on the a- algebra B(Q) 
of the Borel sets of such that for all E G B(Q) the number fi n {E) > yields an estimate, 
in macroscopic average terms, of the amount of people contained in E at time n > 0. Notice 
that this is not a pointwise information on the distribution of pedestrians in ^, but rather a hint 
toward their localization within the walking area, which makes this setup conceptually meaningful 
in terms of a macroscopic look at the system. 

Given the distribution of pedestrians in ft at time n, i.e., the measure ji n : B(Q) — > R + , the 
dynamics of the system toward the next time n + 1 is described by a motion mapping j n : ft — > £7, 
i.e., a Borel function of the form 

j n (x) = x + v n (x)At, (3.1) 

where v n : ft — > R d is the velocity field of pedestrians and At > is the time step. In practice, 
7 n (x) G fl is the position occupied at time n + 1 by the point which at the previous time n was 
located in x G Vt. The new distribution of pedestrians at time n + 1 is obtained pushing the 
measure \i n forward by means of the motion mapping 7 n , which corresponds to defining the new 
measure /i n +i : B(Q) — > as 

V n+1 (E)= f i n ( 1 - 1 (E)), VEeB(n), (3.2) 

or, equivalently written, /i n +i = 7 n #Mn- This equation states in formal mathematical terms the 
simple idea that the amount of people located in a spatial region E at time n + 1 is related to the 
analogous amount at the previous time n along the trajectories of the motion of the pedestrians 



themselves. This interpretation also points out the conservation law structure of Eq. (3.2), which 
can be made even more evident by rewriting the latter in the equivalent form: 

fi n+1 (E) - fi n (E) = -K( 7 - 1 (^ c ) HE)- Hn{l-\E) n E c )\. (3.3) 

Considering that 

1 -\E c )C\E = {x&E : ln (x) #E}, ^(E) n E c = {x $ E : ln (x) € E}, 



we see that the time variation of the measure of E, expressed by the left-hand side of Eq. (3.3) 



is related to the difference between the outgoing flux jj Jn ( / y n 1 (E c ) D E) and the incoming flux 
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p n ( / y~ 1 (E) D E c ), which is the same conceptual idea underlying classical conservation laws in 
continuum mechanics. 

In continuum mechanics one often deals with densities of some continuous fields, thanks to 
which one can switch formally from the integral to the pointwise, namely differential, form of 
balance equations. In our measure-theoretical setting we can recover the concept of density of 
pedestrians by assuming that at time n the measure p n is absolutely continuous with respect 
to the Lebesgue measure C d {continuum hypothesis), written p n <C C d . Then Radon-Nikodym 
theorem implies the existence of a function p n G p n > a.e. in such that d\i n = p n dx, 

namely 

Pn(E) = J p n (x) dx, VE e B(Sl). 

E 

Notice that the pointwise values of p n in Q are physically meaningless, as pedestrians are actually 
not a continuous matter. What instead makes sense is the information on the localization of 
pedestrians and on the crowding of the walking area provided by the measure p n . However, 
having a density is useful not only from the conceptual point of view, due to the analogy with 
classical approaches in continuum mechanics, but also from the numerical point of view, since it is 
much easier to look for numerical approximations of the function p n rather than of the mapping 
jjL n . The question then arises how to guarantee that, given p n <C £ d , also /i n +i is absolutely 
continuous with respect to C d . The answer is furnished by the following 

Theorem 1. Let a constant C > exist such that 

£ d ( 7 - 1 (^)) < C£ d (E), VE e (3.4) 

If ^ £ d then also /i n +i <C C d . In addition: 

(i) ||Pn+i||i = llPnlli, arid 

(ii) if p n e L 1 ^) D L°°(Q) then p n+l e D L°°(Q) as well, with ||p n +i||oo < CWpJ^. 



For the proof of this result the interested reader is referred to [28 . Property (3.4) amounts 
to requiring that j n does not map Lebesgue-nonnegligible subsets of Q into Lebesgue-negligible 
subsets of in other words, that it does not cluster the measure p n during the push-forward. 



It can be proved (see again [28] for the detailed calculations) that j n complies with Eq. (3.4) if 
the velocity field v n is Lipschitz continuous in Q with Lipschitz constant L < At -1 , which can be 
expressively written as: 

At\v n (y) - v n (x)\ <\y - x\, Vx, y eVt. 

If all of the motion mappings {7 n }n>o fulfill this requirement, and moreover an initial measure 
po C d is prescribed, then we deduce from Theorem [l] that p n <C C d all n > with 

llPnlll = lUlll, IKIloo < C"|U||oo. 

It is worth pointing out that, owing to the above-cited Radon-Nikodym theorem, each density p n 
is unique in L 1 ^), whence the well-posedness of Problem ( |3.2[ ) follows, along with the previous a 
priori estimates, for a given initial density po £ ^(Q) D L°°(0). 

Remark (Outgoing flows). Sometimes in applications the motion mappings 7 n may not comply 
strictly with the requirement 7 n (fi) C ^, because for instance the velocity field v n points outward 
the domain in correspondence of some outlet portions of the boundary. For pedestrian flows, 
this commonly happens when exits are present along 917, as we will see more specifically in Sect. 
[5] Nevertheless, the theory illustrated in this section still works, up to some minor technical 
modifications, provided one confines the attention to the restriction measures p n iSl: 

(p n ^n)(E) := p n (E H SI), ME e B(R d ), 

which essentially amounts to neglecting the dynamics possibly taking place outside In partic- 
ular, if 7 n satisfies C d (j~ 1 (E) D Q) < CC d (E) for all Borel sets E C R d , then one gets again the 
existence and uniqueness of densities p n G L 1 ^) D L°°(£}) for the measures _u n \-fl with respect to 
C d . However, the conservation of mass ||p n ||i = ||po||i claimed by Theorem 111 no longer holds, and 
has to be replaced by the more appropriate condition ||p n ||i < ||Po||i- 
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3.2. Numerical treatment of the equations. In order to address the spatial discretization of 
Eq. (3.2) toward the design of a suitable numerical scheme for the approximate solution of the 
problem, we assume conveniently that ft is the cube [0, l] d and partition it using a family of nested 
pairwise disjoint grids {Efyff^l made of cubes and parameterized by the edge size h > of 
the latter: 

n = U E i> Int E i n Int E j = ' 3- ( 3 - 5 ) 

i=l 

In case several, say m, obstacles {Ok}™ = n Ok C SI, are present within the walking area, we further 
suppose that each of them can be described as a suitable union of grid elements, i.e., that for each 
k = 1, . . . , m there exists a collection of indices 1% C {1, . . . , M^} such that Ok = U ieI hE^. In 
principle, the construction described below can be repeated by dropping any, possibly all (except 
properties (3.5)), of these simplifying assumptions, up to additional technicalities in the practical 
implementation of the resulting numerical scheme. 

Let us introduce the measure AJJ <C C d defined by the piecewise constant density 

where \E h * s the indicator function of the set E% , so that dX% = P£ dx. Assume in addition that 
an approximation g 7 ^ of the motion mapping j n is known, having the following structure: 

x\ being the center of the grid cell E^. In essence, this corresponds to an approximation of the 
velocity v n by the piecewise constant function 

M h 
i=l 

which makes a translation over each cell E^ , thus globally a piecewise translation over SI. 

At the successive time step n + 1 we look for a new measure A^ +1 which is in turn absolutely 
continuous with respect to C d and piecewise constant on the partition of Sl\ 

dX n + l = pn + l ^ pn + l W = ^P^Xe^)- 

Imposing that AJJ, A^ +1 satisfy Eq. ( 3.2 ) over each of the grid cells under the action of the mapping 
9h, i.e., 

\ n h +1 (E?) = (g n h #K)(Eh, Vi = 1, . . • , M h , (3.6) 
we obtain, after some standard manipulations, the following scheme: 

Pit 1 = H i = 1, . . . , M h and n > 0, (3.7) 

i=i 

which relates the coefficients of P^ 1 to those of P£ in a time-explicit way. Few observations are 
in order: 

• It can be easily checked that the above scheme is positivity preserving, in the sense that 
if P% > in ft then P£ > in SI all n > 0, and conservative, that is \\P%\\i = ||P£||i all 
n > 0. Moreover, under condition (3.8) discussed below, which, roughly speaking, states 
that the magnitude of the grid size h must be consistent with that of the time step At, 
it is also boundedness preserving: If the initial density is discretized in such a way that 
1 1 -Pj? I |oo < Bo for a certain constant Bo > independent of h, then there exists c > such 
that ||P£||oo < c n B all n > and all h > 0. 
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The image of Ej under g 7 ^ is straightforwardly obtained by exploiting the piecewise trans- 
lation structure of the latter as g^(E^) = E j 1 + v n (x l -)At. 



The measure AJ^ +1 is not the push forward of AJJ by g]\^ indeed Eq. (3.6) holds, by 
construction, on the grid cells only but in principle not for any E G B(Q). It is convenient 
to emphasize this by introducing the notation g\\$- to indicate such a sort of 'approximate 
push forward' operated by the numerical scheme, which however coincides with the actual 
push forward when tested on the grid cells. In other words, g%$- is an operation on a 
piecewise constant measure which returns as output another piecewise constant measure 
(on the same grid as the input measure), whose density with respect to C d is identified 



by a (3.6)-like formula. According to this notation, the scheme rewrites compactly as 

K +1 =9 n JK, n>0. 
The accuracy of the approximation provided by this numerical scheme is related to the ful- 
fillment of a CFL-like condition between the time step At and the grid size ft, as stated by the 
following: 

Theorem 2. Let j n be invertible, smooth, with smooth inverse, and let At, ft > satisfy 

At max \v n (x%)\<h, (3.8) 

i=l, ...,M h 

then: 

(i) (One-step stability) There exist constants A, B > 0, independent of h, such that 

iMn+iW - K + \^)\ < A\\ Pn - PZ\\i + Bh. 

2 = 1 

(ii) (Multistep stability) For each n > there exists a constant A n > 0, independent of h, such 
that 

. max |/in(-Ef ) - \%(Ei)\ < max | Mo (^) - \° h (E?)\ + A n h d . 

i=l,...,M h i=l, ...,M h 

This stability result, proved in detail in [28], allows to control 

(i) the total variation of the localization error in one time step, and 

(ii) the maximum localization error after n time steps 

produced by the approximate measures AJJ with respect to the actual measures \i n in terms of the 
quality 

(i) of the approximation either at the previous time step or on the initial datum, and 

(ii) of the spatial discretization (grid size ft). 

We observe that hypothesis (3.8) of Theorem [2] amounts to controlling the maximum displacement 



of a grid cell, say E^, produced by g\\, indeed for x G we find \g\l(x) —x\ = At\v n (x^)\ < ft. As 
a consequence, the number of nonempty intersections {E% D 9h(Ej)}j%. m Eq. ( |3.7[ ), say c > 0, 
is fixed and independent of the grid size ft (hence also of the total number of cells Mh). Rather, it 
is determined uniquely by the dimension d of the geometrical space in such a way that only grid 
cells Ej sharing an edge with E^ : including obviously E% itself, are at most involved. The number 
c is precisely the constant appearing in the estimate of the oo-norm of P% mentioned above. 

In the case d = 2 this number is 9. Using for convenience a double index notation {E^ k } for the 
grid cells and all the related quantities, this means that only nine adjacent pairs of indices (j, I) 



contribute at most to the intersection with the grid cell (i, k) in Eq. (3.7), namely 

(j, = (i, fc), (2 db 1, fc), (i, k ± 1), (i ± 1, k ± 1), 

so that, denoting by U±, U2 the horizontal and vertical components of v n (x I J l ), respectively, the 
coefficients C 2 (E? k D g^E^)) can be duly computed as 

C\E^ n<(^)) = [U+AtSj^ + (h- p^Sji + C/fA% i+1 ] 

x [E#Attf,, fc _i + {h- \U 2 \At)5 lk + ^ 2 -At<5,, fc+ i], 
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where 5 rs is Kronecker's delta and (-) + ? (*) denote the positive and negative part of their argu- 
ments. 

4. A MODEL FOR PEDESTRIAN MOTION 

As it is clear from the discussion above, the main modeling task to obtain a description of 



pedestrian flows from the measure-theoretical framework outlined in Subsect. |3.1| is the conception 
of a suitable velocity field v n . In particular, the latter must capture some of the most relevant 
behavioral features of pedestrian motion, which specifically differentiate human crowds from other 
systems possibly describable by the same mathematical tools. 

As pointed out by several Authors (see e.g., Bellomo and Dogbe [3], Buttazzo et al. [4], Colombo 
and Rosini [6]), a preliminary important distinction has to be made on the basis of the emotional 
state of the crowd: One can have a flow of pedestrians in either normal or panic conditions, 
each of these states giving rise to dramatically different behavioral rules that people conform to 
collectively. All Authors agree neither on the same qualitative characterization of panic nor on 
the internal mechanisms to the crowd triggering the transition from normal conditions to panic. 
However, at a sufficiently general level of description, we may say that the main difference lies 
in the kind and the strength of interactions that pedestrians experience with one another when 
entering the panic state. In our model we will not be specifically concerned with panic issues. 
Rather we will discuss in detail the main aspects of crowd motion in normal conditions, showing 
how some basic principles, together with an original modeling approach, allow to catch interesting 
features of the system and provide potentially powerful tools for applications in the field. 

A group of pedestrians in motion within a (bounded) walking area Vt C M? usually has a 
target to reach. This may be an exit, a passage, an aggregation point, that we can identify 
with a portion of the boundary dQ., possibly reduced to a point. Therefore the field v n must 
include a first component, call it desired velocity Vd, describing the preferred direction of motion 
of pedestrians toward their target. When the walking area is scattered with obstacles, e.g. pillars, 
walls, furnishings, forbidden areas, which is the most common and interesting case in applications, 
the desired velocity must also account for that pedestrians may need to bypass some of them in 
order to get to the target. 

The desired velocity gives the trajectory that each pedestrian would follow toward her/his target 
in the absence of other neighboring pedestrians. However, an individual walking within a group 
of other pedestrians is influenced in her/his motion by the presence of the latter, mainly because 
she/he will tend to avoid the discomfort caused by highly congested zones. Specifically, pedestrians 
may agree to deviate locally from their preferred path looking for uncrowded surrounding areas. 
Consequently they may decide to correct their desired velocity by an interaction velocity v n , which 
they determine on the basis of the crowding of the walking area in their vicinity. 

The superposition of the desired velocity and the correction operated by the interaction velocity 
yields finally the actual velocity v n of pedestrians: 

4.1. The desired velocity. The desired velocity Vd : Q — > R d advects each individual in the 
walking area toward her/his target, in the (virtual) condition of absence of other pedestrians, 
taking into account that she/he must bypass intermediate obstacles possibly present along the 
path. As such, this velocity field depends uniquely on the geometry of the domain Q and is not 
affected by the distribution of pedestrians within the latter (see Maury and Venel [27]), Vd = Vd(x). 
Some basic technical assumptions we can formulate on the field Vd are the following: 

Assumption 1. Let the desired velocity field v d = Vd(x) : ft — > R d be: 

(i) Lipschitz continuous in Q: 

3T d >0 : T d \v d (y)-v d (x)\<\y-x\, Vz, y G Jl. 

(ii) Uniformly bounded away from zero in Ft: 

3V d >0 : \v d (x)\ >V d , Wxen. 
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The assumption of Lipschitz continuity is related to the homogeneity of the movement of pedes- 
trians in normal conditions. The constant T d , whose inverse is the Lipschitz constant of can 
be interpreted as a characteristic walking time of pedestrians. On the other hand, the uniform 
boundedness away from zero rules out the presence of steady desired flow in Pedestrians always 
have a preferential direction of motion to take. 

If Q is star-shaped with respect to a point xq representing the target of pedestrians, so that 
every x G ^ is connected to xo by a straight path, then a very simple form of the desired velocity 
is obtained as (see e.g., Coscia and Canavesio [8]) 

x - x 

v d {x) = a- 



■ x \ 

where a > denotes the (possibly dimensionless) characteristic magnitude of Vd- We observe that 
such a desired velocity field complies with Assumption [I] up to an arbitrarily small neighborhood 
of xo G Q. 

In more complex domains the above construction of v& is no longer applicable, especially when 
all points xGO are not connected to the target xo by a direct path due to the presence of obstacles. 
Then one may consider the construction of the desired velocity proposed by Maury and Venel [27 , 
which entails the identification of some intermediate targets (often coinciding with some corners of 
the obstacles) to be preliminarily reached until the final target xo is directly accessible. However, 
such a field Vd is in general not as smooth as suggested by Assumption [I] above, because of possible 
sharp swerves when approaching an obstacle. In addition, its generation gets soon laborious for 
an increasing number of obstacles or for variations in their mutual dispositions. 

The method we want to pursue here is based on the introduction of a scalar potential u : Q — > R, 
which identifies attractive and repulsive zones for pedestrians. We assume that u satisfies the 
Laplace equation: 

Au = in tt (4.1) 
and recover the desired velocity field of pedestrians from its (normalized) gradient as 

^7u( x ) 

v d (x) = a a>0. (4.2) 

\Vu(x)\ 

It will be occasionally useful in the sequel to refer to the unit vector z>d(%) — Vu(x)/\Vu(x)\ to 
identify the direction of the desired velocity. 



Equation (4.1) needs to be supplemented by proper boundary conditions in order to yield a 
unique solution u. In our applications we will consider essentially the two most common types of 
boundary conditions, namely: 

(i) Dirichlet boundary condition, consisting in prescribing the value of the potential u along 
the boundary of the domain. 

(ii) Neumann boundary condition, which assigns instead the value of the normal derivative 
|^ along the boundary, n being the outward normal unit vector to dVt. We observe that 



|^ = Vu • n, therefore, in view of Eq. (|4.2|), this kind of boundary condition amounts to 



controlling, in a way, the normal component of Vd on dQ. 

According to the type of boundary under consideration, different conditions might in principle 
apply to different portions of dCl. The ultimate choice often depends on the specific geometry of 
the domain and on the application at hand, however, at a sufficiently general level of description, 
we can state the following guidelines: 

• The outer boundary dft out of ft usually defines the walls delimiting the walking area, 
which may contain doors or, more in general, exits that pedestrians aim at. The latter 
form a subset of dVt ou t that we denote by Ft- To translate the will of pedestrians to reach 
these zones while keeping away from walls, one sets the potential u to zero along the walls 
and to a positive value in correspondence of the targets, for instance: 

(u = ondn out \T T 
1^ = 1 on Tt- 
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Notice that, in order to define a desired velocity (4.2) leading pedestrians to the exits, the 
specific potential of IV is irrelevant provided it is positive, for the gradient \7u will point 
then toward Tt due simply to the growth of u in that direction. 

The inner boundary dQi n is possibly determined by the edges of the obstacles scattered 
within the walking area, that pedestrians must bypass when trying to get to T T . In 
principle, one can model their repulsive effect by prescribing a Dirichlet boundary condition 
of the form 

u = on dQ in . (4.4) 



Coupled with condition (4.3) above, this gives rise globally to a desired velocity field 



pointing outward the obstacles and the outer walls and inward the targets. This is a 



simple consequence of the maximum principle for the Laplace equation (4.1 ), which, owing 
to such a set of boundary conditions, states < u < 1 a.e. in ^, so that the direction of 
Vu is from dVt ou t \TtU dQi n , where the potential attains its minimum, to Tt, where the 
potential attains instead its maximum. 



It is worth noticing that the boundary condition (4.4) usually produces a highly repul- 
sive effect from the obstacles (see Fig. [I] left), while in certain situations, and depending 
on the kind of obstacle, pedestrians may agree to flow by them, simply avoiding head- 
on collisions. This can be rendered by controlling the normal component of the desired 
velocity: 

Vu • n = on dVt in (4.5) 

while allowing free tangential sliding of pedestrians along the walls of the obstacles (Fig. 
[TJ right). 



We remark that the potential u generated by Eq. (4.1) is smooth for domains with smooth 
boundary (and we can assume this is our case, up to possibly round off slightly the sharp corners 
of £1 and of the obstacles), hence Vu is in particular Lipschitz continuous in Q. In addition, u 
being a harmonic function, it results |Viz(x)| 7^ for a.e. x G IntO, for u can have in Int^ at 
most saddle points, due to topological reasons when imposing Dirichlet boundary conditions (cf. 
Fig. [TJ left). It follows that Vd is in turn Lipschitz continuous almost everywhere in Int ^, thus 
finally our desired velocity Vd complies almost everywhere with the requirements of Assumption [I] 
Notice that in case of Neumann boundary conditions, and at least for relatively simple geometries 
of the walking area, we do not expect instead saddle points either (cf. Fig. [I] right). In any case, 
saddle points feature instable manifolds, which prevents the appearance of spurious aggregation 
points within the walking area. 
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v d (x) 



v d (x) 



B+(x) 




\[Mn](x) 



Figure 2. Left and center: The interaction neighborhood B^{x) for a maxi- 
mum angle of visibility max = 7r/2 and < 6 max < 7r/2, respectively. Right: 
Construction of the velocity v n [fi n ](x) as the superposition of the desired velocity 
Vd(x) and the interaction velocity v n \p n ](x) . The point x* is the center of mass 
of pedestrians within B^(x). 

4.2. Interaction velocity. The interaction velocity v n : Vt — > R d expresses the fact that pedes- 
trians are disposed to review and modify locally their preferred trajectory looking for uncrowded 
surrounding areas, because they feel uncomfortable when too close to one another. Unlike the 
desired velocity, v n does not only depend on the geometry of the domain but also on the current 
distribution of pedestrians in i.e., on the measure /i n , which we emphasize by using the notation 

= v n [n n ). Specifically, we assume that each pedestrian evaluates the occupancy of a certain 
neighborhood of her/his position x G ^, call it the interaction (or visibility) neighborhood, and 
then corrects the desired direction Vd(x) by steering clear of the most congested zones she/he finds 
within that neie: hborhood (Fig. |2j. 

Let Br(x) denote the d-dimensional ball of radius R > centered at x G then we take the 
interaction neighborhood of pedestrians located in x G Q to be (Fig. |2| left) 



that is the half-ball in the direction of Vd(x). This corresponds to the idea that pedestrians behave 
mainly like anisotropic particles, in the sense that they look ahead but not behind within a certain 
distance R from themselves. By allowing r(x, y) • Vd(x) > cosO max in the definition above, where 
f(x, y) is the unit vector in the direction of y — x and max G (0, 7r/2], one gets as interaction 
neighborhood the angular sector of B^(x) depicted in Fig. Reenter, with max representing the 
maximum visibility angle of pedestrians (cf. Banos and Charpentier [2 , Bellomo and Dogbe [3]). 

We model the tendency of pedestrians to avoid congested areas by first detecting the center of 
mass x* of the neighborhood B^(x): 



B+(x) 

and then assuming that the correction v n [ii n \ to Vd amounts to a removal from interpreted as 
an indicator of the location of average maximum crowding (Fig. [21 right). We set then: 



where p u [fi n ] : Q, — > is related to the strength of the interaction among pedestrians. Some 
possible choices of this mapping are: 

(i) p u [/^n]( x ) — P/^n(Bft(x))/R for (3 > (conveniently understood as dimensionless), which 
corresponds to a mass-dependent interaction (the more the people the stronger the inter- 
action). 

(ii) p u [fi n ] = P/R for again (3 > 0, which gives a mass-independent interaction (few people 
count as much as many people, the maximum strength of the interaction being fixed by 
/?, indeed in this case |^ n [Mn](^)| < ft all x Ed). 



bU x ) = {y^ b r( x ) • (y-x)- vd(x) > o}, 





(4.6) 
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Figure 3. Correction of a desired velocity field v& (left) by the superposition of 
the interaction velocity v n \p n \ to give the final velocity v n (right). 

Using the regularity hypotheses on Vd expressed by Assumption [I] it is possible to prove that 
option (i) yields a Lipschitz continuous interaction velocity field v n \p n \ (the same is not true 
for option (ii), as Lipschitz continuity of v n [Hn) is broken by the presence at the denominator 
of jj Jn (B^(x)) J which can be arbitrarily close to zero for particular choices of the distribution of 
pedestrians). In short, one first proves that for x±, x 2 G and under the inductive assumption 
that n n is bounded from above by C d (which is, for instance, the case p n <C C d with p n G 
L 1 ^) D L°°(Q) discussed in Sect. |3|, the estimate of \^ n [^n]( x 2) — v n \jjL n ](xi)\ can be essentially 
reduced to the estimate of C d (B^(xi)AB^(x 2 )). Then one shows that for xi, X2 sufficiently close 
the vectors Vd(xi), Vd(x 2 ) point in like directions, that is Vd{x\) • Vd(x 2 ) > 0, which, roughly 
speaking, means that the neighborhoods B^(xi), B^>(x 2 ) almost coincide. This allows to control 
the Lebesgue measure of the non-overlapping parts in terms of the distance between x\ and x 2 as 

C d (B+{ Xl )AB+(x 2 )) <K\x 2 - Xl \ 

for a constant K > independent of xi, x 2l which finally yields the desired Lipschitz estimate on 
v n [n n ]. Analogous calculations are developed in detail in [28] for an interaction velocity built on 
the simpler interaction neighborhood Br(x). 



Notice that the interaction velocity (4.6) gives rise to a nonlinear closure for the velocity v n of 
pedestrians, which in turn depends on \i n (cf. Fig. |3|: 

V n \p>n](x) = V d (x) + V n \fl n ](x). (4.7) 

If, for a given n > 0, fi n <C C d and v n [p n ) is Lipschitz continuous, then we can apply the 
theory developed in [28] for nonlinear fluxes fi n v n [l^n] and recover inductively the well-posedness 
of Problem ( |3.2| ) expressed by Theorem [l] in terms of existence and uniqueness of the densities 
p n G L 1 ^) D L°°(f2), along with a priori estimates on their L 1 - and L°°-norms. 



Remark (Correction to v n near the boundaries of the domain). Equation (4.6) implicitly considers 
in the integration only the portion of B^(x) contained in i.e., B^(x) D ^, because supp p n C ft 
entails that obstacles and perimeter walls are /i n -negligible sets. However, it can be questioned 
that \i n = on ^ c leads erroneously to treat obstacles and walls as uncongested areas, giving 
then rise to an interaction term v n [n n ](x) pointing strongly in their direction, especially for points 
x G Q sufficiently close to dfl (namely, such that B^(x) D ^ c ^ 0). In some situations, this may 
counter the effect of the desired velocity Vd(x) to such an extent that, when summing to obtain 
Vn[l^n](x), one gets a velocity field pointing outside the domain. To correct this drawback it is 
sufficient to compute the interaction velocity with respect to a new measure fi n such that: 

df£ n = 

which extends \i n outside ft by prescribing an equivalent (or effective) density to obstacles and 
walls. The parameter M > can be tuned on the basis of the discomfort produced on the crowd 
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Figure 4. Desired velocity field and initial condition for the problem of the 
narrow passage (i) and of the two narrow passages (ii). 



by the vicinity to the walls, and may possibly vary from wall to wall. Notice that fi n is absolutely 
continuous with respect to C d whenever \i n is, with density p n = p n Xti + Mxn c - 
The expression of v n modifies then in 

vM(x)= P f^\ I (x-y)dUy), (4.8) 

Mn(#fl(aO) J 

B+(x) 

which, from the modeling point of view, renders the fact that obstacles and walls contribute to 
the discomfort felt by pedestrians, in such a way that people might prefer to stay close to one 
another rather than being compressed against the walls. 

Remark (Retrograde flows). The interaction velocity v n [Hn] is, by construction, opposite to the 
desired velocity Vd- It may therefore happen that the resultant v n [ii n ] is in turn opposite to Vd, 
which gives rise to pedestrians locally (in space and time) walking away their targets (retrograde 
flow). If the choice p u [p n ] = P/R is made, it is possible to rule this possibility out by controlling 
the relative magnitude of the parameters a, /3 under the requirement v n [p n ](x) • z>d(x) > all 
xgD, which implies: 

ol + v n [ii n ]{x) ■ v d (x) > 0. 

Since Cauchy-Schwartz inequality implies \v n [p n ](x) • i)d(x)\ < |z/ n [/i n ](x)| < /?, we have that the 
above relation certainly holds true if the following more restrictive condition is satisfied: 

a-0 > 0, 

whence the simple criterion stating that the magnitude of Vd must be greater than the maximum 
possible of f n [^n}' 



5. Application to some cases study 



In this section we apply the time-evolving measures model to some cases study, in order to 
test its ability in reproducing complex realistic features of two-dimensional pedestrian flows. To 



summarize, the model consists of Eq. (3.2) set in Q = [0, l] 2 , supplemented by an initial condition 
yielding the distribution /jlq of pedestrians at time n = 0. The velocity field v n is of the form (|4.7|). 
In particular, we deduce the desired velocity v d from Eq. (|4.2 



solving the Laplace equation (4.1) 
for the potential u along with suitable boundary conditions that will be discussed from time to 



Pfi n (B+(x))/R. 



time, and the interaction velocity v n [{i n ] from Eq. (4.8) choosing p v [fi n ](x) 

In Subsects. 4.1, 4.2 we have discussed the Lipschitz continuity of both Vd and ^ n [Mn]> which, 
together with the assumption /x <C £ 2 , ultimately enables us to apply Theorem [l] and speak of 
density of pedestrians p n G L 1 (^) D L°°(Q) each n > 0. Therefore, in the next applications we 
will constantly refer to these densities when showing the results of the numerical simulations. 
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(i) (ii) (iii) (iv) 

Figure 6. Evolution of the density of pedestrians at successive time instants for 
the problem of two narrow passages. 

5.1. Pedestrian flow through a narrow passage. In this first application we describe the 
motion of a crowd through a bottleneck formed by two long and thin obstacles, parallel and close 
to one another. The bottleneck may represent, for instance, a corridor or some stairs driving 
people toward an exit IV located at its end. The domain of the problem, along with the related 
desired velocity field of pedestrians, and the initial distribution of the crowd are shown in Fig. 
|4^i). In particular, for the computation of Vd we have set u = on all outer boundaries of the 
walking area and on the internal boundaries of the bottleneck, so as to get a repulsive effect that 
leads pedestrians to walk away from perimeter walls and to occupy the middle of the bottleneck. 
Conversely, we have prescribed Vu • n = on the remaining walls of the bottleneck to render a 
minor repulsion by the edges confining with its entrance, and u = 1 on Tt to identify the target 
of the crowd. 

In view of the initial condition depicted in Fig. |4ji), this problem can be compared with an 
analogous application at the microscopic scale proposed by Maury and Venel [27] , who model the 
flow of three groups of pedestrians directed toward a narrow escalator of a station after exiting a 
train. 

Figure [5|i) shows that the three groups join in one single big group few instants after the 
beginning of the simulation, and then walk toward the entrance of the bottleneck from the left. 
Once there, they give rise to an obstruction (Fig. |5jh)), with some people pushed also on the right 
side of the entrance (Fig. |5jiii)), until they all flow through the bottleneck (Fig. [5|iv)). Notice 
how the model predicts a self-organization of pedestrians in lanes of high density separated by 
areas of lower crowding. This is especially evident in Fig. [5|ii) for the group of pedestrians sliding 
along the upper edge of the first obstacle, and also in the bifurcation of the flow occurring in front 
of the entrance of the bottleneck. 

5.2. Pedestrian flow through two adjacent narrow passages. We address now a case similar 
to that dealt with in Subsect. |5.1[ but we add a second bottleneck beside the previous one, that 
equally leads to an exit (Fig. |4Jii)), so that pedestrians have the choice of which passage to 
take in order to reach their target. Starting from the initial condition depicted in Fig. (4^ii), the 
model predicts an evolution initially similar to that of the previous problem (compare Fig. |5|i) 
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Figure 7. Density of pedestrians in a domain with two obstacles, starting from 
the initial condition (i). The desired velocity field is that depicted in Fig. [I] with 
(ii) Dirichlet and (iii) Neumann boundary conditions for the potential u at the 
edges of the obstacles. Snapshots (ii) and (iii) refer to the same evolution time 
instant. 



and Fig. [6^i)). Later on, some individuals, pushed sideways by the crowd at the entrance of the 
first bottleneck, decide to take the other passage (Fig. (6^ii)), which gives rise to a flow also in 
the second bottleneck. Pedestrians then start branching off in correspondence of the separation 
between the two bottlenecks (Fig. |6jiii)), until all of them have entered either passage (Fig. [6^iv)). 

This problem may model, for instance, the flow of passengers exiting the cars of a subway and 
then heading for the exit of the metro station. Many stations have two (or even more) adjacent 
escalators that pedestrians can use to reach the ground level, however it is commonly observed 
that, on average, they prefer to take the closest one to their starting point. This behavior is clearly 
caught by the model, indeed Figs. [6jii-iv) show a lower crowd density in the farthest passage. 

5.3. Effect of boundary conditions and modeling of obstacles. Boundary conditions for 
the desired velocity field Vd at the obstacle edges may sensibly affect the resulting configuration 
of the flow of pedestrians. Let us consider a walking area Q featuring two obstacles Oi, O2 and 
an exit T T located on a portion of the right vertical boundary, as illustrated in Fig. [I] A compact 
group of pedestrians is initially positioned on the opposite side of the room (Fig. |7ji)) and is 
guided toward the exit by either the field Vd depicted in Fig. [I] left, with the potential u set to 
zero at every internal and external boundary of Q (but obviously IV, where condition u = 1 is 



ll right, characterized by zero normal derivative 
iii) show an instantaneous configuration of the 



prescribed instead), or the field Vd depicted in Fig. 
of the potential u at the obstacles. Figures [T^ii) , [7 
distribution of pedestrians in either case at the same evolution time. The repulsive effect of the 
obstacles is particularly evident in the case of Dirichlet boundary conditions (Fig. [T^ii) ) , with 
preferential paths bypassing both obstacles from the outside. Neumann conditions (Fig. (7^iii)) 
produce instead the expected sliding of pedestrians along the edges, with a larger crowding of the 
area comprised between the two obstacles. This demonstrates that the set of boundary conditions 
to generate the desired velocity field has in general nontrivial consequences on the resulting flow 
of pedestrians, and may vary from case to case also in connection with the role played by each 
obstacle in every specific application. 

5.4. Lane formation versus clustering. We focus now on the self-organization properties pre- 
dicted by our model for a system of agents characterized by either anisotropic or isotropic inter- 
action mechanisms. Figure |8ji) shows a compact group of agents initially located near the left 
side of the domain Vt. We assume that they are driven by a constant rightward desired velocity 



field, which originates from the potential u(x, y) = x solving Eq. (4.1) with the following set of 
boundary conditions: 

u = on x = 

1 on x = 1 

Vu • n = on y = 0, 1. 
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(i) (ii-b) (iii-b) 

Figure 8. Self-organization of a group of agents. Starting from the same initial 
condition (i), anisotropic interactions lead to lane formation (ii-a, iii-a), whereas 
isotropic interactions promotes clustering (ii-b, iii-b). 

Assume that the agents under consideration are pedestrians, who, as already mentioned, feature 
anisotropic interactions between each other for they are able to see ahead only. As shown in 
Fig. |8jii-a, iii-a), this fosters a self-organization of the mass in parallel lanes, which follow the 
direction of the main velocity field and whose number depends on the initial distribution fio of 
pedestrians. We notice that the reciprocal distance between the lanes is comparable to the size of 
the neighborhood of interaction (R = 0.1 in this simulation). Organization in uniformly walking 
lanes is actually observed in real pedestrian flows, as pointed out by Helbing and coworkers in 

HQunun]. 

Let now the agents be able to see both ahead and behind, as it happens for instance with birds, 
whose view covers also a back area due to the lateral position of their eyes on the head. Interactions 
among the agents are then isotropic, as they do not depend on the desired direction of motion 
but are equally felt on the whole ball Br(x). In such a situation, starting from the same initial 
condition as before, our model predicts a clustering of the mass (Fig. (8^h-b, iii-b)) while drifted 
by the main velocity field, with a characteristic distance between clusters again comparable to the 
size R of the neighborhood of interaction. This kind of cooperative motion pattern is actually 
found in several multi-agent systems, such as birds or fishes, and in the specialized literature is 
commonly termed flocking (cf. e.g., Krause and Ruxton [24] and the main references therein). 

5.5. Crossing flows. Finally, we study the problem of two groups of people walking in opposite 
directions toward one another. Several Authors have addressed this issue, especially from the 
microscopic point of view, finding in general a good agreement between the main features of this 
kind of flow and the predictions provided by their models. In the papers by Helbing and cowork- 
ers [TTJ [14] and by Hoogendoorn and coworkers [T9j HO] ? which are a valuable source of empirical 
information on pedestrian behavior, particular attention is given to typical patterns emerging in 
crossing pedestrian flows. The microscopic models proposed by Helbing et al (ibidem), Hoogen- 
doorn and Bovy [17 , Maury and Venel [26] have proved to be successful in reproducing such 
experimental evidences. 

In order to model crossing flows we need to introduce two measures /i* , /i^ : B(Q) — > R+ 
describing the space occupancy by either pedestrian group, and consequently two motion mappings 
7^, 7^ : Vt — > Vt such that 



^ +1 (E) = (ii((Y n ) \E)), VE e B(fi), i = 1, 2. 



(5.1) 
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(i) 



(ii) 



(iii) 



Figure 9. (i) Two uninterrupted flows of pedestrians traveling rightward (upper 
row, population i = 1) and leftward (lower row, population i = 2) meet at the 
center of the domain. The two groups start to interact with each other, breaking 
the initial left-right symmetry (ii) and giving rise finally to a stable configuration 
of alternate uniformly walking lanes (iii). The cyan line in each plot is the 0.1-level 
curve of the density p % n . 



Each group has its own desired velocity v % d : Q — > R d , which, consistently with the characterization 
stated in Subsect. |4.1[ is unaffected by the desired velocity of the opposite group. Conversely, 
each interaction velocity depends now on both pedestrian distributions: 

K = K^li /4L i = 2 > 

indeed the dynamics triggered by encounters among individual belonging to different populations 
plays a relevant role in this problem. Specifically, we assume that pedestrians of the first population 
feel uncomfortable when too close to pedestrians of the second population, due to their different 
walking targets, and vice versa. Therefore, they may decide to keep away from areas of high 
concentration of people coming in the opposite direction, trying instead to gain room for their 
walking direction. At the same time, people may still disagree with being overcompressed by 
pedestrians of the their own population, and therefore avoid an excessive proximity also with 
those walking in the same direction, or prefer instead to keep the contact with the latter in order 
to predominate over the oppositely walking population. To take these phenomena into account, 
we propose the following form of the interaction velocity: 

1 r 2 

<[^,^](x) = - J (x-y)d^^j^(y), * = 1,2, (5.2) 

B+(x) 3 = 1 



possibly along with an obvious extension of the correction (4.8) at the boundaries of the domain, 
where: 

• Pij > for i 7^ j determines the strength of the interaction of the population i with the 
head-on population j; 

• flu G R gives the strength of the interaction among people of the same population, with a 
disgregating effect if flu > and an aggregating effect if instead flu < 0. 



The mathematical problem generated by the two coupled equations (5.1), supplemented by 



proper initial conditions /ij, /i§ <C £ d , can be profitably addressed by making use of the techniques 



illustrated in [28] for nonlinear flows. In particular, the form (5.2) of the interaction velocity along 



with the regularity of the two (uncoupled) potential desired velocity fields v\, v\ give rise to motion 
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mappings 7^, 7^ both complying with hypothesis (3.4) of Theorem [T] This enables one to prove 
that two sequences of densities 

{/£}»><>, {^WoCL'^nrift) 

exist, each of which is unique, such that p l n > a.e. in ft and dfi^ = p % n dx each i = 1, 2 and each 
n > 0. 

Figure [9] shows the prediction of the model for two uninterrupted flows of pedestrians entering 
the domain from the left and the right boundary, respectively. In this simulation, we have set 
fin — P22 — (no interaction among pedestrians belonging to the same population), so as to 
focus specifically on the effect of the interactions between the two oppositely walking groups. The 
results clearly demonstrate the ability of the model to reproduce first the rupture of the left-right 
symmetry, as soon as the two groups start to interact, and then the emergence of typical patterns 
of alternate uniformly walking lanes, which are experimentally observed as a characteristic self- 
organization phenomenon in crossing flows and, at present, mathematically described mainly by 
means of microscopic models (cf. the above-cited papers). 

It is worth noticing that our mathematical model has not been conceived with the specific 
purpose of describing either the self-organization in lanes/clusters discussed in Subsect. 5.4 or 



the emerging uniformly walking lanes addressed here. Rather, we stress that the model is able 
to catch all these behavioral customs as a by-product of much more general (and even more 
elementary) modeling principles, resorting ultimately to the basic idea of nonlocal interactions 
among pedestrians. 



6. Conclusions and research perspectives 

In this paper we have systematically applied the mathematical structures developed in [28] to 
the modeling of pedestrian flows. The reference framework consists of a discrete-time sequence 
of Radon positive measures {fi n }n>o defined over the Borel a-algebra 23(0), ft C R d being the 
domain of the problem, which evolve according to the recursive relation /i n +i = 7 n #Mn (push- 
forward). The core of the structure is represented by the motion mappings j n : Q — > O, j n (x) = 
x + v n (x)At, which describe the dynamics of the system. Here, v n : Q — > R d is the velocity field 
of pedestrians and At > the time step. The regularity of the measures /i n , specifically the 
fact of being absolutely continuous with respect to the Lebesgue measure C d on R d , depends on 
some non-singularity properties of the 7 n 's, which can be rephrased in that they do not cluster 
the Lebesgue measure, i.e., do not map Lebesgue non-negligible sets into Lebesgue negligible sets. 
This guarantees the existence and uniqueness of densities {p n }n>o C L 1 (0) such that d\i n = p n dx, 
which are also bounded in O, i.e., p n G L°°(Q) each n > 0, provided po is (cf. Theorem [T]). The 
existence of densities is particularly significant for application purposes, as it allows to establish 
a direct connection of the measure-theoretical structures with classical methods of continuum 
mechanics. Indeed, not only can the relation /i n +i = 7 n #Mn be equivalently rewritten in a form 
which reminds of the usual structure of conservation laws (cf. Eq.( |3.3| )), but, considering that the 
push forward of the /i n 's turns out to be the direct time discretization of the abstract continuous- 
time conservation law (cf. [28] ) 

j t »t(lt(E)) = 0, \/Ee B(Q), 

one is led, under the continuum hypothesis /i^, /i n <C C d each t > 0, n > 0, to a time discretization 
of the equation 

j t J p(t,x)dx = 0, VEeB(Q), 

lt(E) 

i.e., the classical (Lagrangian) form of the mass conservation equation dealt with by continuum 
mechanics. 

Numerical approximation of the push forward can be performed by a suitably devised scheme 
(cf. Eq. ( |3.7| )), introducing piecewise constant measures AJJ : B(Q) — > M + over a partition of 
the domain Q with characteristic size individuated by a parameter h > 0. Under some technical 
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assumptions, one can control the so-called localization error (cf. Subsect. 3.2) produced by the 
AJJ's over the /i n 's in terms of the grid size h (cf. Theorem |2j), provided the latter is chosen with 
respect to the time step At in such a way to fulfill a CFL-like condition (cf. Eq. ( |3.8| )). 

It is worth stressing that this theoretical setting enables one to address both analytical and 
numerical issues in a thorough and unified manner for whatever dimension d > 1 of the spatial 
domain. The same is not true for classical approaches based on nonlinear hyperbolic conservation 
laws, which are well known to generate nontrivial analytical and numerical difficulties as soon 
as the spatial dimension of the problem gets greater than 1. This is particularly relevant for 
the application to pedestrian flows, which develop mainly in two dimensions: The time-evolving 
measures framework allows a qualitative mathematical mastery of the problem that nonlinear 
conservation laws might hardly achieve. As an example of the difficulties raised by standard 
modeling techniques, let us consider the handling of boundary conditions. From the modeling 
point of view, one has to guarantee that pedestrians do not either flow through obstacles or leave 
the walking area from portions of the boundary other than the prescribed outlet regions. On 
the other hand, when using hyperbolic partial differential equations, the definition of the outflow 
portion of the boundary is not simply a matter of appropriate boundary conditions, because at 
each time instant one must also take into account the orientation of the characteristic velocities of 
the problem. Therefore, it may be argued that some unilateral constraints on the velocity should 
be imposed at the boundary, that the solution of the problem must then comply with, for instance 
v • n ^ on some portions of dfl, in order to "deflect" characteristics where needed. Clearly, this 
introduces additional technicalities, that can be instead ruled out elegantly by a time-evolving 
measures approach, as indicated in Subsect. |4.2| 



The main modeling task posed by the time-evolving measures setting (3.1 )-(3.2) is the definition 
of the velocity field v n . For pedestrian flows, inspired by some analogous considerations proposed 
by Maury and Venel [27 , we have suggested a structure accounting for two main factors affecting 
the motion of walkers, namely: 

• Pedestrians' will to reach a specific target placed in the walking area, expressed by a 
desired velocity field : ft — > R d , which essentially depends on the geometry of the 
domain, including the presence of intermediate obstacles to be bypassed. This is the 
velocity a person would have in the absence of other people in the surroundings, and is 
modeled as a conservative vector field coming from a scalar potential. 

• Pedestrians' tendency to avoid congested areas, expressed by an interaction velocity v n [/jL n ] 
Q — > R d , which depends in a nonlocal way on the distribution fi n of pedestrians them- 
selves. To model this term, we have borrowed some ideas from the Eulerian approach to 
the analysis of rendez-vous problems for multi-agent systems by Canuto et al. [5]. 

The superposition of these two effects, plus some suitable corrections at the boundaries of Q and 



of the obstacles illustrated in Subsect. |4.2[ yields the final form of pedestrians' velocity: 

V n \p>n](x) = V d (x) + V n [{l n ](x), X G ^, 

which closes the model. 

By means of these structures, we have addressed some representative cases study with the aim 
of testing the ability of the model to reproduce complex features of pedestrian flows pointed out 
in the specialized experimental literature (see e.g., Helbing et al. [TOj HH H3] , Hoogendoorn et al. 
[I9j|20], Krause and Ruxton [24], and their main references). In more detail, we have considered 
several applications to motion in domains with obstacles, possibly also with more than one target 
for pedestrians, mimicking realistic situations like the access to one or two escalators by a group 



of people getting off the cars of a subway (cf. Subsects. 5.1, 5.2) or the outflow of a crowd from a 
room with two pillars partially hiding the exit (cf. Subsect. 5.3). In the latter case, we have also 
studied the effect of different boundary conditions for the desired velocity field at the edges of the 
obstacles on the flow. In addition, we have shown that the model is able to reproduce interesting 
self-organization phenomena of human crowds, like the formation of dynamic lanes in a group of 



pedestrians sharing the same desired direction of motion (cf. Subsect. 5.4) and the emergence of 



alternate uniformly walking lanes in crossing flows (cf. Subsect. 5.5). 
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It may be questioned that, unlike standard macroscopic models of continuum mechanics, for 
instance those using the formalism of partial differential equations, the modeling structures by 
time-evolving measures presented in [28] and in the present paper do not allow to deal with 
continuous-time dynamics. In other words, we are currently lacking a "limit model" for At — > 
0. It is worth pointing out that, as far as applications are concerned, this is actually not a 
dramatic drawback, because the framework we have introduced still provides conceptual tools 
for a mathematical study of pedestrian flow problems, which have proved to be successful in 
addressing numerous applications to real cases. On the other hand, the great reward for the price 
paid of discrete time modeling is the possibility to approach in a unified manner, from both the 
analytical and the numerical point of view, d-dimensional systems with no additional technical 
difficulties when d > 1. We stress the importance of this both to attain a satisfactory theoretical 
mastery of the modeling structures and to treat immediately realistic two-dimensional applications, 
without the need for devising preliminarily simplified explorative one-dimensional approximations. 
In general, the same is not as straightforward for more classical modeling techniques based on 
nonlinear conservation laws. Finally, we observe that the investigation of the limit behavior of the 
model for At — > is, however, a possible by-product of this research line, which has then also the 
merit of promoting new ways of theoretical speculation. 

Besides theoretical issues, further research developments may involve the use of the model of 
pedestrian flows presented here in connection with control and optimization topics. It is known 
that obstacles, which in the present work have been regarded as essentially passive elements of the 
walking area, often play instead an active role in forcing the direction of motion of the walkers. As 
an example, they may be designed and positioned so as to foster the flow of the crowd along specific 
directions, as it happens for instance in shopping centers when buyers are induced to move close to 
particularly attractive zones, or to improve the safety of pedestrians when accessing/leaving certain 
areas. In this context, we would like to recall Braess 7 paradox for pedestrian flows, which has been 
pointed out by Hughes [23] and then drawn by several other Authors. Such a paradox, credited 
to the mathematician Dietrich Braess, was originally formulated for traffic flows on networks: It 
states that adding extra capacity (for instance, a further edge) to a network can in some cases 
reduce the overall performance (see e.g., [31] for more details). Rephrased in the abstract, it means 
that a condition intuitively expected to lead to an improvement may instead give rise to worse 
outcomes. Hughes suggests to invert Braess' paradox for pedestrian flows, after observing that 
placing an obstacle in front of an exit (intuitively worse condition) may sometimes improve the 
stream of people from a room (better outcome) . The reason for this unexpected fact might be due 
to a sort of reorganization of the flow of pedestrians induced by the obstacle, provided shape, size, 
and positioning of the latter are accurately studied. Clearly, the mathematical treatment of such 
issues needs to take advantage of sufficiently handy, but also realistic, two-dimensional models of 
pedestrian flows. In this respect, we believe that our macroscopic model by time-evolving measures 
can be profitably used to ground specific control and optimization problems of practical interest 
for real word applications. 

Other research perspectives can be developed on the modeling side. As reviewed in Sect. [2] 
some macroscopic models already available in the literature (cf. Bellomo and Dogbe [3 and 
the main references therein) describe the flow of pedestrians by both the conservation of mass 
and the balance of linear momentum, using Navier-Stokes-like partial differential equations with 
nonclassical force contributions. A natural question is whether one can do the same in a measure- 
theoretical setting, specifically what could be a measure-theoretical counterpart of the balance 
law 

d t (pv) + V • (pv i;) = F[p, v], 
which, as it is well-known in continuum mechanics, is equivalent to the second equation of system 



(2.2) under the assumption of conservation of the mass. The question gets even more interesting 
if one considers that the nonlinearity of the convective term V • (pv v) makes Navier-Stokes 
equations able to catch the possible appearance of turbulent regimes in the fluid, which Helbing et 
al. [12] pointed out to play a role also in the flow of pedestrians. It may be argued that a careful 
translation of the continuum linear momentum balance equation in the language of time-evolving 
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measures could enable our model to describe also turbulence phenomena in the motion of human 
crowds. The transition from laminar to turbulent flow could then be interpreted as the ability of 
the model to reproduce the emergence of panic. 

References 

[1] A. Aw and M. Rascle. Resurrection of "second order" models of traffic flow. SIAM J. Appl. Math., 60(3) :916- 
938, 2000. 

[2] A. Banos and A. Charpentier. Simulating pedestrian behavior in subway stations with agents. In F. Amblard, 
editor, Proceedings of The Fourth Conference of the European Social Simulation Association (ESSA2007), 
pages 611-621. September 2007. 

[3] N. Bellomo and C. Dogbe. On the modelling crowd dynamics from scaling to hyperbolic macroscopic models. 
Math. Models Methods Appl. Sci., 18(suppl.):1317-1345, 2008. 

[4] G. Buttazzo, C. Jimenez, and E. Oudet. An optimization problem for mass transportation with congested 
dynamics. Preprint, Department of Mathematics, University of Pisa, Italy (http://cvgmt.sns.it), 2007. 

[5] C. Canuto, F. Fagnani, and P. Tilli. A Eulerian approach to the analysis of rendez-vous algorithms. In Pro- 
ceedings of the 17th IFAC World Congress (IFAC'08), pages 9039-9044. Seoul, Korea, July 2008. 

[6] R. M. Colombo and M. D. Rosini. Pedestrian flows and non-classical shocks. Math. Methods Appl. Sci., 
28(13):1553-1567, 2005. 

[7] R. M. Colombo and M. D. Rosini. Existence of nonclassical solutions in a pedestrian flow model. Preprint, 
2008. 

[8] V. Coscia and C. Canavesio. First-order macroscopic modelling of human crowd dynamics. Math. Models 

Methods Appl. Sci., 18(suppl.):1217-1247, 2008. 
[9] P. Dallard, T. Fitzpatrick, A. Flint, S. Le Bourva, A. Low, R.M. Ridsdill, and M. Willford. The London 

Millennium Footbridge. The Structural Engineer, 79(22):17-33, 2001. 
[10] D. Helbing, I. J. Farkas, P. Molnar, and T. Vicsek. Simulation of pedestrian crowds in normal and evacuation 

situations. In M. Schreckenberg and S. D. Sharma, editors, Pedestrian and Evacuation Dynamics, pages 21-58. 

Springer, Berlin, 2002. 

[11] D. Helbing and A. Johansson. Quantitative agent-based modeling of human interactions in space and time. 

In F. Amblard, editor, Proceedings of The Fourth Conference of the European Social Simulation Association 

(ESSA2007), pages 623-637. September 2007. 
[12] D. Helbing, A. Johansson, and H. Z. Al-Abideen. Dynamics of crowd disasters: An empirical study. Phys. Rev. 

E, 75(4):046109, 2007. 

[13] D. Helbing and P. Molnar. Social force model for pedestrian dynamics. Phys. Rev. E, 51(5):4282-4286, May 
1995. 

[14] D. Helbing, P. Molnar, I. J. Farkas, and K. Bolay. Self-organizing pedestrian movement. Environment and 

Planning B: Planning and Design, 28(3):361-383, 2001. 
[15] S. Hoogendoorn and W. Daamen. Self-organization in walker experiments. In Traffic and Granular Flow '04- 

2004. 

[16] S. P. Hoogendoorn and P. H. L. Bovy. State-of-the-art of vehicular traffic flow modelling. J. Syst. Cont. Eng., 
215(4):283-303, 2001. 

[17] S. P. Hoogendoorn and P. H. L. Bovy. Simulation of pedestrian flows by optimal control and differential games. 

Optim. Control Appl. Meth., 24:153-172, 2003. 
[18] S. P. Hoogendoorn and P. H. L. Bovy. Pedestrian travel behavior modeling. Networks and Spatial Economics, 

5:193-216, 2005. 

[19] S. P. Hoogendoorn and W. Daamen. Self-organization in pedestrian flow. In Traffic and Granular Flow '03, 

pages 373-382. Springer, Berlin Heidelberg, 2005. 
[20] S. P. Hoogendoorn, W. Daamen, and P. H. L. Bovy. Extracting microscopic pedestrian characteristics from 

video data. In Transportation Research Board annual meeting 2003, pages 1-15. National Academy Press, 

Washington DC, 2003. 

[21] R. L. Hughes. The flow of large crowds of pedestrians. Mathematics and Computer in Simulations, 53:367-370, 
2000. 

[22] R. L. Hughes. A continuum theory for the flow of pedestrians. Transportation Research Part B, 36:507-535, 
2002. 

[23] R. L. Hughes. The flow of human crowds. Annu. Rev. Fluid Mech., 35:169-182, 2003. 
[24] J. Krause and CD. Ruxton. Living in Groups. Oxford University Press, Oxford, 2002. 

[25] M. J. Lighthill and G. B. Whitham. On kinematic waves. II. A theory of traffic flow on long crowded roads. 

Proc. Roy. Soc. London. Ser. A., 229:317-345, 1955. 
[26] B. Maury and J. Venel. Handling of contacts in crowd motion simulations. In Traffic and Granular Flow '07. 

Springer. To appear. 

[27] B. Maury and J. Venel. Un modele de mouvements de foule. ESAIM: Proceedings, 18:143-152, 2007. 

[28] B. Piccoli and A. Tosin. Time-evolving measures and macroscopic modeling of pedestrian flow. Preprint, 

arXiv:0811.3383vl [math-ph]. Submitted. 
[29] P. I. Richards. Shock waves on the highway. Operations Res., 4:42-51, 1956. 



PEDESTRIAN FLOWS IN BOUNDED DOMAINS WITH OBSTACLES 



25 



[30] F. Venuti, L. Bruno, and N. Bellomo. Crowd dynamics on a moving platform: mathematical modelling and 

application to lively footbridges. Math. Comput. Modelling, 45(3-4) :252-269, 2007. 
[31] Wikipedia. Braess's paradox — Wikipedia, the free encyclopedia, 2008. 

ISTITUTO PER LE APPLICAZIONI DEL CALCOLO "MAURO PlCONE" , CONSIGLIO NAZIONALE DELLE RlCERCHE, VlALE 
DEL POLICLINICO 137, 00161 ROMA, ITALY 

E-mail address: b.piccoli@iac.cnr.it 

ISTITUTO PER LE APPLICAZIONI DEL CALCOLO "MAURO PlCONE" , CONSIGLIO NAZIONALE DELLE RlCERCHE, VlALE 
DEL POLICLINICO 137, 00161 ROMA, ITALY 

E-mail address: a.tosin@iac . cnr . it 



